Ab initio zone-center phonons in LiTaO,3: comparison to LiNbOs 
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The four Ai-TO V phonon frequencies in lithium tantalate are calculated in the frozen-phonon 
approach from first principles using the full-potential linearized augmented plane wave method. A 
good agreement with the experimental data available is found for all modes; reliable displacement 
pattern of different modes becomes available from the calculated eigenvectors. The Raman spec- 
tra recorded for Ai modes in LiNbC>3 exhibit a counter-intuitive softening of the j4i-T03 mode 
frequency with respect to that measured in LiTa03. We explain this behavior by a comparatively 
harder oxygen rotation in LiTa03 and discuss other differences in lattice dynamics of two materials, 
notably derealization of Ta and Li contributions over more that one corresponding mode in LiTa03 , 
differently from the situation in lithium niobate. The Li isotope shift is predicted in the calculation. 
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INTRODUCTION 

The ferroelectric materials LiNbC>3 and LiTaC>3 have 
been intensively studied over years due to their promis- 
ing applications in non-linear optical and electro-optic 
devices - see Ref. @ for a review on relevant properties 
of these materials. Both compounds possess rhombohe- 
dral space group R3c (C® v ) and have 10 atoms in the unit 
cell. They remain ferroelectric up to quite high transition 
temperatures: T c amounts to 1480 K in LNLO3 (that is 
probably the highest ferroelectric transition temperature 
so far known) and 950 K in LiTaOa. 

Due to low symmetry and large number of atoms in 
the unit cell, the first-principles studies of LiNbOs and 
LiTaC>3 proceed not so fast as, e.g., for perovskite-type 
ferroelectric compounds, KNbC>3 and KTaC>3 (the latter 
is an incipient ferroelectric, where the transition can be 
induced by uniaxial pressure, or by dopingl. Soon after 
first ab initio studies of electronic structurea, the analy- 
sis of th£ aspects of ferroelectric instability by Inbar and 
Cohenou was illuminating in that they emphasized an 
importance of oxygen-rotating mode in the ferroelectric 
transition. Earlier, the ^-displacement of Li ions out of 
the oxygen planes was believed to be the essence of the 
ferroelectric transformation. Inbar and Cohen demon- 
strated that the energy profile in the para- to ferroelec- 
tric transition (studied along the path connecting exper- 
imentally determined end-point structures) is very sim- 
ilar in both substances. As a possible explanation of 
the difference in the T c values, they guessed that the 
zone-boundary behavior may be found to be different in 
LiNbC>3 and LiTa03. This assumption has not yet been 
tested. 

While we do not access exactly this topic in the present 
study, we believe that the analysis of phonon properties, 
even at zone center, can give a clue to understanding 
different behavior of LiNbC>3 and LiTaOs. Due to the 
complexity of structure, the analysis of lattice vibrations 



in experiment was in part contradictory, and in what 
regards theory only few calculations have been done. 
For LiNbOs we calculated earliertffl the frequencies and 
eigenvectors of TO-r phonons, based on a precision total- 
energy fit in the multidimensional (up to dimension 9, for 
tha.E block) space of symmetry coordinates. Parlinski et 
alxi calculated the phonon dispersion over the full Bril- 
louin zone of LiNbOs, using the direct method and the 
Fourier transformation of force constants sampled by dis- 
placing individual atoms in a 80-atoms (2x2x2) super- 
cell. Semiempirical phonon calculations (i.e., fitting the 
model parameters to experimentally measured frequen- 
cies) have been done by Repelin et aln for both LiNbOs 
and LiTaOs. We are not aware of any first-principles 
phonon calculations for LiTaOs. 

On the experimental side, at least for the A\ modes 
there is very good agreement between different measure- 
ments in what regards vibration frequencies - see Table @. 
Comparing them with the frequencies for LiNbOs shown 
in the same table one can note that the frequencies are 
lower in LiTaOs f° r all modes but one, the TO3. Whereas 
the mode softening can be - at least qualitatively - intu- 
itively related to larger mass of Ta, and different degree of 
softening of different modes - to different participation of 
Ta (Nb) in them, the hardening of the TO3 mode on go- 
ing from LiNbOs to LiTaOs is clearly counter-intuitive 
and can be only traced to the differences in electronic 
properties. 

The aim of the present study is to provide an ab initio 
description of zone-center Ai-TO phonons in LiTaOs and 
discuss them in comparison with corresponding phonons 
in LiNb03, thus outlining small but important differ- 
ences in electronic properties in these materials. In the 
following we briefly outline the details of calculation and 
proceed with the discussion of phonon frequencies and 
eigenvectors. 
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TABLE I. Ax-TO phonon frequencies in LiNbQ 3 and LiTa0 3 
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TABLE II. Internal coordinates in LiTa0 3 (experiment and theoretical optimization) 







Ta 






Li 






O 




Experiment: Ref. |17| 
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Present calculation 
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0.049 
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0.071 



CALCULATION SCHEME 

Our frozen-phonon calculation was based on precision 
all-electron total-energy calculations with the use of full- 
potential linearized augmented plane waves method, as 
implemented in the WIEN97 codaij. The local orbitals 
extension according to Ref. |lj has been used for a bet- 
ter description of Ta 5s and 5p, Li Is and 02s states. 
The muffin-tin radii chosen were 1.88 a.u. for Ta and 
1.65 a.u. for Li and O. The exchange-correlation energy 
has been treated in the local density approximation as 
parametrized by Perdew and WangEi 

The Brillouin zone integrations—were carried out by 
the improved tetrahedron methocO on a 4x4x4 special 
fc-points mesh which generated 13 k points in the irre- 
ducible Brillouin zone. When using a denser 6x6x6 k- 
mesh (32 irreducible fc-points), we found a difference in 
the total-energy trends less than 1 mRy, negligible for the 
analysis of the lattice dynamics in the scale of energy vari- 
ations discussed below. The planewave cutoff parameter 
-Rmin^max for the basis was set to 9.5, that resulted in, 
on the average, 2313 basis functions for each fc-point; the 
cutoff parameter G max for the charge density expansion 
was 12.0, i.e., these cutoff values had to be set tcuhigher 
values than in our earlier calculations for LiNbOsQ in or- 
der to achieve acceptable convergence. The optimized 
lattice parameters turn out to be a=5.1274 a.u. and 
c=13.7033 a.u., giving the error of —2% in the volume 
and —0.1% in the c/a ratio as compared to experimental 
values a=5.1543, c=13.7835 a.uO. The subsequent cal- 
culations analyzing the atomic displacements have been 
performed with the lattice parameters fixed at their ex- 
perimental values. 

As regards the iipzen-phonon formalism as such, it is 
known since longjiS and has been successfully applied 
to many systems. This approach is physically trans- 
parent, and its accuracy is easy to control. The bot- 
tleneck of the method is the necessity to calculate the 
forces (or, alternatively, expansion coefficients of the 
total energy over individual atomic displacements) for 
many trial geometries, in order to map the total energy 
near equilibrium to quadratic or, if necessary, higher or- 



ders in displacements. This rapidly makes the compu- 
tation unpractical for moderately large systems, and/or 
those with reduced symmetry (where many vibrational 
degrees of freedom mix in the same irreducible represen- 
tation). Earlier, we succeeded in calculating zone-center 
phonons far, at the largest, the 9-dimensional S-block 
of LiNbOaQ, adding progressively data points (i.e., trial 
geometries) that would refine normal coordinates, as the 
latter emerge more and more clearly in the course of cal- 
culations. In the present study, our concern is the 4- 
dimensional A\ block of the crystal vibrations in LiTaOa, 
for which we make a comparison with previously obtained 
LiNb03 dataB. In order to obtain a reliable mapping of 
the total energy we included 170 different geometries in 
our calculation. The optimized internal coordinates as 
they emerge from the four-dimensional total-energy fit 
are shown in Table H They serve as reference point for 
all displacements discussed later on. 



RESULTS 

For the set of four symmetry coordinates S1...4 chosen 
in terms of individual Cartesian displacements as speci- 
fied in Ref. || for LiNb03, the diagonalization of the dy- 
namical matrix yields four normal coordinates as follows: 

f S TO! \ 
S T0 2 
S TO s 
V S T0 4 ) 

These normal coordinates serve as abscissa in Fig. [j] 
that visualizes corresponding cuts of the total energy hy- 
persurface. The number of trial configurations allows to 
perform a stable 4th order polynomial fit, also shown in 
Fig. [|. It is obvious that the corrections to the 2d or- 
der fit are small and not always physically sane, i.e. not 
all leading 4th order terms are positive. In other words, 
all four Ax-TO modes, according to our calculation, are 
essentially harmonic, much more so than it was the case 
for their counterparts in LiNbOs (see a similar analy- 
sis in Ref. 0). Another "empirical" justification of this 
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FIG. 1. The cut of the total-energy fit along the directions of normal displacements in the four-dimensional space of 
symmetry coordinates corresponding to the Ai symmetry. The circles' sizes are inversely proportional to the distances of 
corresponding data points from the direction of normal displacement. Solid lines represent the 2d-order polynomial fit, dashed 
lines - the 4th order fit. 



fact is much better agreement of measured and calculated 
(harmonic) frequencies than it was the case in LiNb03 
(see Table |). The largest difference (7 - 12 cm -1 ) be- 
tween our calculation and experiment occurs for the TOi 
and TO2 modes, with experimental values higher, as it 
could be usually expected due to anharmonicity. For two 
other modes, the agreement between experiment and the- 
ory is almost perfect, as consistent with neat parabolic 
fit through data points along the directions of these two 
normal coordinates in Fig. |l|. In LiNb03, the difference 
between theory and experiment amounted to 46 cm -1 
(TOi) and 50 cm -1 (TO4 mode), both these modes being 
considerably unharmonic judging on their corresponding 
total energy profiles. In the following, we discuss the 
composition of different modes and address the "anoma- 
lous" TO3 mode (i.e., that is harder in tantalate than in 



niobate) in some detail. 

Since the LiNb03 and LiTaOs are in many aspects so 
similar, we would like, when discussing their phonon fre- 
quencies, to know whether the existing differences are 
mostly due to higher tantalum mass, or to variations in 
corresponding force constants. The simplest illustration 
of the effect of mass is the calculation of frequencies with 
the mass of Nb atom taken instead of that of Ta; the 
results are given in the last column of Table |[ Unex- 
pectedly enough, the frequency noticeably increases in 
the course of such substitution solely for the TO2 mode; 
for three other modes, even for the TOi (the mainly 
"Nb" mode in LiNb03, see Ref. ||) the effect is negli- 
gible, consistently with much smaller Ta contribution in 
corresponding eigenvectors (see below). 

Such "Ta isotope" effect roughly accounts for the half 
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Li-z displacement (a.u.) 

FIG. 2. Potential energy profile associated with Li 
z-displacements in A1-TO3 mode of LiNbOs (LNO) and 
LiTaC>3 (LTO). Positive movements shorten the distance be- 
tween Li and O ions, while the negative movements shorten 
that between Li and Nb(Ta) ions. 

of the TO2 frequency shift in Li(Ta— >Nb)C>3; the rest 
can only come from the steeper potential well associated 
with Li z-displacement in LiNbC>3. The corresponding 
total energy profile is shown in Fig. ||. The reason for 
the "softness" of the LiTaC>3 lattice with respect to such 
Li-only vibration can be simply a somehow larger inter- 
atomic spacing, as illustrated by Table ffl. It was already 
emphasized earlier by Inbar and Cohenff that Li always 
remains quite well ionized in both tantalate and niobate 
and that no effects of dynamical covalency that would 
lead to coupling of oxygen and lithium motions are seen 
in these materials. 

Let us now turn to calculated eigenvectors shown in 
Table |v| in the hexagonal setting. The "prototype" oxy- 
gen atom is that whose coordinates are given in Table [fi|; 
from its eigenvector components (a:, y, z) the components 
related to other oxygen atoms can be reconstructed as 
(y,x-y,z), (y-x,x,z), [y-x,y,\ + z), (y,x,\ + z), 
(x, x—y, i + z). For comparison, the eigenvector compo- 



TABLE III. Distances (a.u.) to nearest neighbors and next 
nearest neighbors of transition metal and oxygen to lithium 
atom in LiNb03 and LiTaC>3 





Nb(Ta) NN 


Nb(Ta) NNN 


Onn 


Onnn 


LiNbOs 


5.595 


7.159 


3.799 


4.276 


LiTaOa 


5.669 


7.355 


3.848 


4.367 



nents for LiNbC>3 are also given; their transformation to 
the Cartesian system according to 
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recovers the eigenvectors for LiNbC>3 as given in in Ref. |6|. 
As in LiNbC>3, the heavy transition- metal ion contributes 
mostly to the two low-frequency modes. In the softest 
one, it vibrates along with Li against the oxygen octa- 
hedra which are only slightly tilted in this process (the 
tilting originates merely from a crystal structure frustra- 
tion, as was discussed in Ref. ||). This is the ferroelec- 
tric mode, frozen down from the paraelectric phase, near 
its "freezing point" . The TO2 mode includes lithium 
vibration against both transition-metal and oxygen sub- 
lattices; two other modes are too high for efficiently res- 
onating with transition- metal vibrations. 

The Li contribution exhibits some differences in two 
systems. To make these differences, somehow obscured 
in center-of-mass related eigenvectors, more clear, we 
show in Table |y| pure (not scaled with masses) vibra- 
tion patterns of the four modes, relative to the heaviest 
Nb(Ta) atom and normalized so that the displacement of 
(any) second heaviest constituent atom, i.e. oxygen, is 1. 
Whereas in LiNbC>3 the lithium vibration is essentially 
confined to the^TX^ mode (that alone is affected by the 
Li isotope shiftBH), in LiTa0 3 the Li effect on the TOi 
and TO3 modes is also not negligible, with possible im- 
plications for their frequencies on Li isotope doping. The 
frequencies calculated with reduced Li mass are shown 
in Table |. The maximal frequency shift (for the TO2 
mode) reduces from 19 cm -1 in LiNbC>3 to 10 cm -1 in 
LiTaC>3. We guess that the Li isotope shift of the TOi 
mode is probably still large enough to be detected ex- 
perimentally. Moreover, the Li-related phonon density 
of states, if ever calculated similarly to how it was done 
in Ref. for LiNbOs, must certainly become spread off 
lower and higher frequencies. 

The second largest contribution of Li occurs, for both 
materials, in the eigenvector of the TOi mode. Even if 
it is much smaller (relative to the center of masses) in 
niobate, the displacement pattern of this mode is oxygen 
z-vibration against both the heaviest atom and lithium, 
so that the Li-0 stretching is anyway present. Here again 
the hardening of the Li-z-related potential well on chang- 
ing from LiTaC>3 to LiNbC>3 plays a role: the mere change 
of mass produces only a minor effect in elevating the cal- 
culated (harmonic) frequency to 208 cm^ 1 in niobateEl; 
the strong anharmonicity of this mode contributes the 
rest up to its measured frequency of 252 cm -1 . 

The two hardest modes exhibit but negligible Nb(Ta) 
contribution in both materials; the effect of Li z- 
displacement is also quite reduced (although remains no- 
ticeable in LiTaOs), so we have here two almost exclu- 
sively O xy patterns, almost identical in tantalate and nio- 
bate. These two modes are, correspondingly, a nearly 
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TABLE IV. Eigenvectors of TO phonons in LiNbOa and LiTaOs in the hexagonal setting 



LiTaOii LiNbOa 





Ta 


Li 




O 




Nb 


Li 




O 




TOi 


0.060 


0.168 


( 0.096 


0.132 


-0.104 ) 


0.143 


0.035 


( 0.036 


0.085 


-0.123 ) 


T0 2 


0.109 


-0.191 


(-0.047 


0.038 


-0.080 ) 


0.068 


-0.254 


( 0.015 


-0.014 


0.001 ) 


TO3 


-0.007 


0.051 


(-0.441 


-0.101 


-0.003 ) 


0.007 


-0.001 


(-0.463 


-0.156 


-0.006 ) 


TO4 


0.029 


0.022 


(-0.127 


-0.440 


-0.037 ) 


0.023 


0.015 


(-0.082 


-0.437 


-0.022 ) 



TABLE V. Cartesian coordinates at equilibrium (top line) and displacement patterns in the four Ai modes of LiNb03 and 



LiTaC>3, relative to the Ta(Nb) atom and with the oxygen displacements normalized to 1. 





Li 




01 






02 






03 
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z 
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y 
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LiTaOs 












coord. 


0.282c 


-0.280a 


0.122a 


0.405c 


0.034a 


-0.304a 


0.405c 


0.246a 


0.181a 


0.405c 


TOi 


0.75 


0.13 


-0.04 


-0.38 


-0.04 


0.13 


-0.38 


-0.10 


-0.10 


-0.38 


TO2 


-1.14 


0.04 


0.09 


-0.40 


-0.10 


-0.01 


-0.40 


0.05 


-0.08 


-0.40 


TO3 


0.22 


-0.09 


0.40 


-0.00 


-0.30 


-0.28 


-0.00 


0.39 


-0.12 


-0.00 


TO4 


0.07 


-0.38 


-0.09 


-0.12 


0.27 


-0.28 


-0.12 


0.11 


0.37 


-0.12 












LiNbOa 












coord. 


0.281c 


-0.278o 


0.127a 


0.401c 


0.029a 


-0.304a 


0.401c 


0.249a 


0.177a 


0.401c 


TOi 


-0.01 


0.06 


0.01 


-0.40 


-0.04 


0.05 


-0.40 


-0.03 


-0.06 


-0.40 


TO2 


-6.17 


-0.07 


-0.05 


-0.40 


0.07 


-0.03 


-0.40 


-0.01 


0.08 


-0.40 


T0 3 


-0.01 


-0.14 


0.38 


-0.02 


-0.27 


-0.31 


-0.02 


0.40 


-0.08 


-0.02 


TO4 


0.04 


-0.38 


-0.14 


-0.08 


0.31 


-0.26 


-0.08 


0.07 


0.39 


-0.08 



rigid rotation of oxygen octahedra hosting Nb atoms and 
the bending of the octahedra (rotation of three top atoms 
against the three bottom ones) . The "top view" of cor- 
responding distortions has been shown in Fig. 2 of Ref. || 
for LiNb0 3 . What is noteworthy about T0 3 and TO4 
modes is that for both of them the calculated frequency is 
by ~16 cm -1 higher in tantalate. Then, the anharmonic- 
ity of TO4 in LiNbOa (but not in LiTaOs) accounts for 
a considerable increase of experimental frequency in the 
niobate. The TO3 mode is well harmonic in both ma- 
terials, so that the relation between phonon frequencies 
in niobate and tantalate as predicted by the calculation 
is essentially retained in reality. Since the TO4 mode 
in both materials is almost free of transition-atom and 
lithium contributions, the last panel of Fig. fil and the 
corresponding energy profile in Fig. 1 of Ref. [3| exhibit 
essentially the effect of bending of oxygen octahedra and 
give a fair indication that the latter are more rigid in 
LiTa0 3 . 

An analogy to that can be found in cubic per- 
ovskites: among many calculations done by date of zone- 
center phtppps in KNbOa and KTa03 by comparable 
techniqueaiaEd, the calculated frequencies of two hardest 
modes, that correspond to different distortions of oxygen 
octahedra, are higher in tantalate than in niobate. The 
experimental situation confirms this at least for the TO4 
mode (the TO3 of the T^ u symmetry is Raman and in- 
frared silent in cubic perovskite, so that its measurements 
are not numerous; corresponding frequencies any anyway 
quite close in KNb0 3 and KTaOs). The possible origin 



of this seemingly systematic trend can be sought for ei- 
ther in a slightly increased lattice spacing of tantalates 
- that would however probably have an opposite effect 
on the oxygen rigidity, if at all - or in a slight difference 
in the degree of covalency. The latter guess gets support 
from the analysis of the dynamical matrix corresponding 
to the soft-mode motion (i.e., also including the oxygen 
tilting) by Inbar and CohciJJ. They have shown that the 
"full" 0-0 elements of this matrix are larger in LiTa03, 
whereas the relation between Madelung contributions to 
the dynamical matrix (i.e., those assuming nominal ionic 
charges) is opposite. 

The aspects of relative covalency of tantalates and nio- 
bates are influenced by several factors; discussions to the 
point can be found in Ref. ^ on perovskites and in Ref. |] 
on our present systems of interest. Summarizing, the co- 
valency seems to be higher in tantalates, judging by a 
slightly wider valence band and a stronger dispersion of 
certain bands, that is seemingly related to less localized 
d-function of tantalum. A possible link to the stiffness of 
oxygen octahedra may be that higher covalency enhances 
interaction between lattice-dynamical degrees of freedom, 
softening the softest modes and hardening the hardest. 
Indeed, the span of TO-frequencies, at least in the har- 
monic approximation, seems to be typically broader in 
tantalates. Different degree of anharmonicity of different 
modes may complicate the experimental situation. 
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CONCLUSION 

We extended our previous frozen-phonon analysis of 
zone center-TO phonons in lithium niobate over lithium 
tantalate. For all A\ modes quite good agreement was 
found between the experimental and calculated phonon 
frequencies, indicating, differently from LiNbC>3, a quite 
harmonic nature of all four modes. In spite of large 
similarity in crystal structure and electronic proper- 
tics, the phonon eigenvectors show certain differences 
in both systems. Specifically, the Li contribution, that 
was in LiNb03 only noticeable in the single TO2 mode, 
now spreads over TOi to TO3 modes. One can ex- 
pect the Li isotope shift in LiTaC>3 to become less pro- 
nounced but affecting these three modes. Li displace- 
ment experiences a softer potential well in LiTa0 3 , due 
to more expanded lattice, that lowers the frequencies of 
Li-dependent modes. The contribution of the transition 
metal ion affects the TO2 mode more than TOi, con- 
trary to the case in LiNbC>3. Contrary to the effect on 
Li, the rotational modes of oxygen octahedra are notice- 
ably hardened in LiTa03, indicating a more stiff oxygen 
network. Due to this, a higher frequency of the TO3 
mode in tantalate is explained. For the TO4 mode, the 
harmonic part of which also gets harder in LiTaC>3, the 
experimental situation is reversed only due to a strong 
anharmonicity of this mode in LiNbC>3. 
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